Anisotropies in the Gravitational Wave Background from Preheating 
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We investigate the anisotropies in the gravitational wave (GW) background produced at preheat- 
ing after inflation. Using lattice field theory simulations of a massless preheating model, we show 
that the GW amplitude depends sensitively on the value of the decay product field \ coupled to the 
inflaton (j>, with the only requisite that \ is light during inflation. We find a strong anisotropy in 
the amplitude of the GW background on large angular scales, the details of which strongly depend 
on the reheating dynamics. We expect similar conclusions for a wide class of inflationary models 
with light scalar fields. If future direct detection GW experiments are capable of detecting the GW 
produced by preheating, they should also be able to detect this effect. This could eventually provide 
a powerful way to distinguish between different inflationary and preheating scenarios. 



The last two decades have seen enormous progress 
in cosmology, driven mainly by measurements of the 
temperature anisotropies of the cosmic microwave back- 
ground (CMB) radiation, and culminating in the recent 
publication of the data from the Planck satellite |l| . In 
principle, gravitational waves (GW) have the potential to 
become an equally important source of information about 
the early universe. Their great advantage is that while 
the universe was opaque to electromagnetic radiation for 
the first ~ 400.000 years, GW were able to travel freely 
through space ever since their emission, and therefore 
they can carry information to us from the earliest mo- 
ments of the universe [2J|, close to the Big Bang. Possible 
sources of GW in the early universe include quantum fluc- 
tuations during inflation |3| , non-equilibrium phenomena 
after inflation [J, |fj| , and cosmic defects [6| . 

One particularly interesting source of GW is reheating, 
the transition from inflation to a radiation phase. During 
this process, the inflaton field, which was driving the in- 
flation, decays into other fields. In many models, reheat- 
ing begins with an initial non-equilibrium stage known as 
preheating pj , which is much more sensitive to the details 
of the model than the period of inflation itself. In par- 
ticular, preheating generates a sub-horizon background 
of GW with a significant amplitude |5| , whose frequency 
depends on the energy scale of inflation, being typically 
in the MHz-GHz range. This is too high to be observed 
with existing or planned large GW detectors such as Ad- 
vanced LIGO, but prototype detectors for 100 MHz GW 
have been built [8|- 

In the existing literature, the focus has been on de- 
termining the spectrum of the GW, but the purpose of 
this letter is to show that under certain relatively general 
conditions, the GW amplitude can depend on position on 
cosmological scales. This means that the observed GW 
background will be anisotropic on large angular scales. 
Such anisotropy should be potentially observable with 
any GW experiment that is capable of detecting the GW 
background, and might provide a powerful new way to 



constrain models of inflation and preheating. 

The anisotropy arises in models with a second light 
scalar field x, in addition to the inflaton field. By light we 
mean that its mass is less than the Hubble rate H* dur- 
ing inflation, thanks to which it should develop an almost 
Gaussian scale-invariant spectrum of fluctuations during 
inflation (similarly to the inflaton), with power spectrum 
V x sw H%/4tt 2 Q. When inflation ends, the details of 
preheating in each Hubble volume depend on the local 
value of the \ field, which we denote by Xi- I n particu- 
lar, this will also be true for the energy density p GW of the 
GW produced within each Hubble volume, which there- 
fore becomes a local function of Xi- That is, normalizing 
the GW energy density to the critical energy density p c 



today, fl c 



,j Pc, we expect O c 



^Gw(Xi)" 



The GW background from preheating originated from 
a comoving spherical shell of radius R ~ l/i?0j with Hq 
the Hubble rate today. This surface encompasses a very 
large number of preheating Hubble patches. The GW 
background we could detect in a given direction n, was 
then produced within one of these patches centered at 
r = Rn. We expect therefore that its amplitude should 
depend on the value of Xi within the preheating Hub- 
ble volume at r = Rn. Because Xi nas fluctuations on 
cosmological scales, so does fi GW , and therefore the ob- 
served GW amplitude will depend on the direction n, 



i.e. 



n GW (h) = n GW (xi(Rn)). 



The details of the anisotropy are determined by how 
the GW energy density depends on xu i- e - by the function 
!1 GW (x0- We will present a lattice calculation of f2 GW (x0 
for a massless preheating model shortly, but to illustrate 
the effect let us first assume a linear dependence as 



^Gw(Xi) = CQ + Cl((5Xi/-ff*), 



(1) 



where cq and c\ are dimensionless constants, and Sxi = 
Xi — Xi i with X; the average value over the currently ob- 
servable universe. In this case the relative fluctuations of 
the GW energy density are <5fi GW = (fl GW /fl GW - 1) = 
(ci/co)(Sxi/H*)- Like 5xi, these fluctuations are nearly 



Gaussian and scale-invariant, and their power spectrum 



is V r _ 



( Cl /c o y(P x /Hi) = (ci/27rc ) 2 . The angular 



power spectrum of the relative fluctuations on large scales 
(i.e. small multipole I), is given by 
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(2) 



in analogy with the Sachs- Wolfe plateau [9| . 

In practice, ^ Gw (xi) wu l not be linear and therefore 
the calculation of the effect is more complicated [10j • In 
order to quantify the anisotropics on large angular scales, 
we compute the two-point correlation function of the GW 
energy density originated at two points x and y. The 
joint probability distribution for the field values Xx = 
Xi(x) and Xy = Xi(y) a t these points is 



■P(Xx,Xy) «exp 



<^l(Sxl + $xl) - 2G^ y 5 X x6xy 



2(4 



^x,y) 



(3) 

where Sx = Xi ~ Xi> G x, y = ( S Xi( x )$Xi(y)) is the field 
correlator, and ct 2 = G xx = (Sx 2 ) is the field variance. 
The GW energy density correlator is then given by 

(tt GW (x)tt GW (y)) = 

dXxdXy-P(Xx,Xy)^ G w(Xx)ft GW (Xy)- (4) 



The analysis simplifies if the correlator is well approx- 
imated by its linear Taylor expansion in powers of the 
field correlator G* 



z x,y, 



(n GW (x)ft GW (y))^(tt GW ) 2 + 



($x& aw (xi)Y 



G x ,y, (5) 



where the expectation values on the right hand side are 
computed with the single-point probability distribution 

P(Xi)cxexp{-( X i-Xi) 2 /24}. (6) 

Then, the GW energy density correlator coincides with 
the linear ansatz ([I]) with coefficients c$ = (r2 GW ) and 
c\ = -ff* (Jxi^Gw(Xi)) /fly' The angular power spectrum 
is therefore given by Eq. @, 

Hi (5 Xi n GW ( Xi )) 2 



l(l + l)Q 



8ir a$(n c 



(7) 



We have calculated this explicitly in the massless pre- 
heating model. It is one of the easiest models to analyze 
numerically, and it has been extensively studied analyti- 
cally and numerically 111, |l2| . It describes the inflaton </>, 
with a chaotic quartic potential during inflation, coupled 
to a scalar \ as 



mx) = ^ 4 + ^W, 



with A = 9 x 1CT 14 as fixed by the amplitude of the CMB 
anisotropics. This model is conformally invariant, which 
makes it very convenient for lattice simulations, since 
comoving physical scales will remain within the lattice 
volume throughout the simulation. Although massless 
preheating has been nearly ruled out by observations [13JJ , 
we use it for its numerical convenience. We expect our 
conclusions to be valid for more realistic models as well. 

During inflation <p is large, <p ^ Mpi, and x is as ~ 
sumed to be very small. When inflation ends, <f> starts 
oscillating around its minimum at (j) = 0. These oscil- 
lations induce an instability in the fluctuations of the 
field x- The amplitude of the x modes initially grows ex- 
ponentially fast for certain comoving momentum bands 
as Xk(t) oc e Mfc, '^- ) , where drj — dt/a(t) is the conformal 
time. This effect is known as parametric resonance [7| 
and has been studied analytically in great detail for this 
model in [11| . The unstable solutions where Xfe grows ex- 
ponentially are characterized by an exponent fik, which 
depends on the wavenumber k and the ratio of coupling 
constants g 2 /X- The band structure of this model is sum- 
marized in Fig. 4 of 11|. We focused on the parameter 
choice g 2 / X — 2, for which x is light during inflation and 
parametric resonance is strongest for the longest wave- 
lengths. Eventually the exponential growth ends when 
the linear approximation breaks down. 

The GW are primarily produced during the subsequent 
stage of non-linear non-equilibrium dynamics, which we 
study by using lattice field theory simulations. These 
consist of solving a discretized version of the scalar field 
equations of motion of the matter fields in the model, 

4> + 3H<j>~ ^V 2 + (\4> 2 + g 2 X 2 )^ = 0, (9) 

X + 3H X - 4 V 'X + 9 2 4> 2 X = 0, (10) 

coupled to the Friedmann equation 3a/a = — 47r(p + 
3p)/3Mpj for the scale factor a(i), sourced by the volume 
averaged pressure p and energy density p contributed by 
(f> and x- The Planck mass is defined in terms of Newton's 
constant G as Mpi = G -1 / 2 . Our code was based on the 
publicly available MPI C/C++ ClusterEasy package fl4J . 
which solves the field equations with a second-order leap- 
frog integrator using periodic boundary conditions. 

To compute the GW spectrum we followed the algo- 
rithm proposed in 15( . During preheating the system de- 



velops an anisotropic stress tensor Ily = T, 
Tij is the energy-momentum tensor. 



The 



pgij, where 
transverse- 



traceless part of Ily acts as a source for the tensor pertur- 
bation hij representing GW. In practice, the GW source 
is given by ITJ T = {dixdjX + di<j)dj(j)} TT , where {• • • } TT 
stands for a projection of the transverse-traceless (TT) 
degrees of freedom, guaranteeing 9iII^ T = II J T = 0. 
Therefore we solve a discretized version of the linearized 
Einstein equations for the tensor perturbations hij in the 
l°J presence of the non-linear field evolution determined by 
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FIG. 1: Final spectrum of GW for Xi = 3.4 x 10" 8 A/ P1 (up- 
per, blue curve) and Xi — 1-0 x 10 _8 Mpi (lower, black curves), 
averaged over five random realizations of inhomogeneous fluc- 
tuations. The solid curves are for L — f60, N — f024, and 
the dashed curves for L — 80, N = 512. The area underneath 
corresponds to the total fractional GW energy density within 
a preheating Hubble domain. 



Eqs. @ and (HUD, 



3Hh t 



'■h lJ = 
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W 



FT 



all^^x). (11) 



We verify the conditions dihij — hu = by choosing a 
lattice-based TT-projector associated to neutral lattice 
derivatives [16(, but checked that the choice of projec- 
tor did not influence our results. Finally, we obtained 
the total energy density of GW within a volume V and 
normalized to the critical energy density, as 



"owW 



1 

f><- 



( d Pc 



dp G 



k 3 M^ 



dlogk (At:) 3 V 



\ d log k 
dflk 



dlogk 



(12) 



-in 



%(£,fc,k)/i?-(*,A:,k) (13) 



For the discretized version of Eq. (fT3"|) . see Eq. (4.5) 
in [16|. The GW amplitude grows with time, and even- 
tually saturates, reaching a final amplitude at some time 
tf, when the fields enter a turbulent regime. 

As explained earlier, we want to determine how the 
final fractional GW energy density f2 GW (£/) generated 
within a single Hubble region of volume V — H~ 3 de- 
pends on the initial value Xi within same the region. 
Therefore, we repeated this calculation for a range of 
different Xi- This is in contrast with most of the liter- 
ature on massless preheating to date, in which Xi has 
been assumed to be zero. The main exception is the cal- 
culation curvature perturbations produced by preheat- 
ing, which depends on the impact of Xi on the expansion 
rate [T7( . These studies have demonstrated a chaotic be- 
haviour which leads to non-Gaussian spikes in the curva- 
ture perturbation C [18|. 



FIG. 2: fi GW for our sample of initial field values Xi- 

The initial conditions for the fields <f> and x consisted 
of inhomogeneous Gaussian fluctuations mimicking vac- 
uum quantum fluctuations [19j |, superimposed on homo- 
geneous background values (j>o = 0.342Mpi, xo = Xi- 

Our simulations show the expected behavior of para- 
metric resonance: the inflaton oscillates initially, trans- 
ferring energy to the fluctuations in x which grow until 
the system becomes non-linear. As the fluctuations grow, 
so does the gradient energy of the fields which leads to 
GW production. The field gradients eventually become 
an adiabatic slowly evolving function, which sets the end 
of production of GW. In practice, the GW spectrum con- 
tinues to oscillate, and we therefore took a time average 
over several oscillations to obtain the final GW ampli- 
tude. Fig. Q] shows the shape of the GW spectrum for 
two initial field values Xi- The dashed^ lines corresponds 
to our fiducial choice of lattice size L = vA<A)£ = 80 
and lattice points per dimension N — 512, whereas the 
solid lines correspond to L = 160 and N — 1024 [ensur- 
ing the same ultraviolet (UV) coverage]. For L = 160, 
one can clearly see a drop in the infrared (IR), which 
shows that very long wavelength modes are not excited. 
In practice, the total integrated GW amplitude for the 
two volumes agree to better than ~ 1%. Because L = 160 
was computationally too expensive and it was not pos- 
sible to capture both IR and UV behaviors sufficiently 
well with smaller lattices, we chose L = 80, iV = 512 for 
our simulations. 

To determine the anisotropy of the GW background, 
we need to compute the expectation values in Eq. ([7]) 
with the Gaussian distribution ([6]), with variance 



H 



r H * dk 

a l = / TT^x = f^^CMB « 3.3 x 10- ib Af^, (14) 

Ja H k 47T 



where we have used the Hubble rate at the end of inflation 
Hi k 87rA$J/12M£i « 2.6 x 10~ 15 M^ and iVcMB ~ 50 
as the number of e- folds after the largest observable scales 
left the horizon. Because it is likely that inflation lasted 



more than 50 e-foldings, the scale invariant fluctuations 
of the x held continue outside our current horizon, and 
therefore we have to allow for a non-zero average value 
Xi over the currently observable universe, with variance 
(Xi) = {H 2 jAir 2 ){N tot - iVcMB), where N tot is the total 
number of e-foldings. We treat Xi as a f ree parameter 
with approximate magnitude Xi ~ 10~ 7 Mpi. 

We calculated the expectation values in Eq. ((7]) using 
the Monte Carlo method, generating JV = 400 random 
values Xi j J £ {1, . . . ,N"} from the Gaussian distribution 
©, choosing the mean value Xi = 3-42 x 10~ 7 Mpi. For 
each Xi, we did one simulation run, measuring the GW 
energy density ^l G w(Xi )' see Fig- 121 As the plot shows, 
f2 GW is highly dependent on xu varying by as much as 
a factor of five between nearby values, although there 
are some ranges of Xi where the dependence is much 
smoother. This irregular behavior is in line with the 
chaotic dynamics observed earlier 18[ , but its amplitude 
is unexpectedly high. 

We computed the full correlator (J4} using the Monte 
Carlo data to confirm the validity of the linear expansion 
([5]) and Eq. ([7]). The expectation values in Eq. are 
approximated by averages within our sample, 



0.03 






if 2-~l GwlXi / 



N 



Tf^2^d -Xi)^Gw(xi'. 



N 



(15) 



For Xi = 3.42 x Kr 7 M P i, we obtained (0 GW ) = (5.45 ± 
0.13) x 10~ 4 and (Sx^ GW ) = (3.0 ± 1.2) x 10- 12 M P1 . 
Substituting these into Eq. (J7]) gives the amplitude of 
the relative fluctuations <5S1 GW = (r2 GW /fi GW — 1) as 



y/l{l + TjCi = 0.017 ± 0.008, 



(16) 



where the errors are estimated by the bootstrap method. 
We extended our results to nearby mean values x! by 
reweighting our data, giving to each x{ the weight 



exp 



(xf-xO 2 , (xl-x) 2 



M 



2ct 2 



(17) 



in the averages (TT5j) . with x and a 2 the actual numer- 
ical average and variance of our sample. The resulting 
amplitudes are shown in Fig. [3] From these data we can 
conclude with some confidence that the relative ampli- 
tude of the anisotropies is above one per cent level in 
this case, much higher than the CMB anisotropies. 

In summary, we have shown that the GW background 
from preheating is generally anisotropic on large angular 
scales if light scalar fields are present during inflation. 
We have used numerical lattice field theory simulations 
to demonstrate that in the massless preheating model 
the amplitude of these anisotropies is significant, at the 
~ 1% level. Obviously, measuring these anisotropies will 




FIG. 3: The relative amplitude of the multipoles of the 
GW background as a function of the average field value Xi, 
calculated from Eq. (0. The black dot shows the amplitude 
for original mean value Xi ~ 3-42 x 10 _7 Mpi, and the curve 
shows values obtained by reweighting the same data. 



be a great challenge even if the GW background itself 
is detected, but if achieved, it might provide detailed 
information about the microscopic physics of inflation. 

The massless preheating model is exceptional because 
of its conformal invariance, and it is likely that in more 
realistic theories the GW energy density f2 GW is less sen- 
sitively dependent on the field value Xi> like the curvature 
perturbation [20(. This does not necessarily mean that 
the amplitude of the anisotropies will be lower, because 
in the current case the highly irregular variation in Fig. [2] 
becomes largely averaged out. Our general conclusions 
are also not necessarily restricted to models with para- 
metric resonance, and we would generally expect the GW 
background to be anisotropic on large angular scales in 
any model in which GW emission is affected by a light 
scalar field. 



In light of the recent discovery of a Higgs-like particle 
at the LHC [21( , it would be interesting to investigate this 
phenomenon in a realistic particle physics scenario, with 
the Standard Model Higgs field playing the role of the 
light scalar field. This will be computationally more de- 
manding and will require the inclusion of gauge fields 22] 
and fermions 12311. 
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